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ABSTRACT 

We present the first study that combines binary population synthesis in the Galactic disk and detailed evolu- 
tionary calculations of low- and intermediate-mass X-ray binaries (L/IMXBs). This approach allows us to follow 
completely the formation of L/IMXBs, and their evolution through the X-ray phase, to the point when they become 
binary millisecond pulsars (BMPs). We show that the formation probability of IMXBs with initial donor masses 
of 1.5— 4 Mq is typically >5 times higher than that of standard LMXBs with initial donor masses of <1.5M Q . 
Since IMXBs evolve to resemble observed LMXBs, we suggest that the majority of the observed systems may 
have descended from IMXBs. Distributions at the current epoch of the orbital periods, donor masses, and mass 
accretion rates of L/IMXBs have been computed, as have orbital-period distributions of BMPs. Several signifi- 
cant discrepancies between the theoretical and observed distributions are discussed. We find that the total number 
of luminous (Lx > 10 36 ergss _1 ) X-ray sources at the current epoch and the period distribution of BMPs are very 
sensitive to the parameters in analytic formula describing the common-envelope phase that precedes the formation 
of the neutron star. The orbital-period distribution of observed BMPs strongly favors cases where the common 
envelope is more easily ejected. However, this leads to a > 100-fold overproduction of the theoretical number of 
luminous X-ray sources relative to the total observed number of LMXBs. X-ray irradiation of the donor star may 
result in a dramatic reduction in the X-ray active lifetime of L/IMXBs, thus possibly resolving the overproduction 
problem, as well as the long-standing BMP/LMXB birthrate problem. 

Subject headings: binaries: close — pulsars: general — stars: neutron — X-rays: stars 



1. INTRODUCTION 

Roughly 140 low-mass X-ray binaries (LMXBs) are known 
in the Galaxy (Liu et al. 2001), with orbital periods from 1 1 min 
to ~1 yr, donor masses of ~0.01-2M Q , and X-ray luminosities 
from the detection sensitivities to ~10 38 ergss _1 . Over the past 
twenty years or so, many theoretical studies of LMXBs have 
aimed at accounting for their abundance and variety. During 
this time, a standard picture for the formation and evolution 
of LMXBs in the Galactic disk has emerged. However, recent 
observational and theoretical work have challenged the conven- 
tional wisdom and prompted a renewed interest in the origins of 
observed LMXBs. Specifically, it has been realized that many, 
perhaps even the majority, of the identified LMXBs with low- 
mass stellar companions may be descendants of systems with 
intermediate-mass (>1.5Mq) donor stars. 

In the past, all binary population synthesis studies that ex- 
plicitly considered the evolution of X-ray binaries and the cri- 
teria for dynamically unstable mass transfer involved analytic 
approximations (e.g., Rappaport et al. 1982; Kalogera & Web- 
bink 1996; King & Ritter 1999). This is a satisfactory ap- 
proach as long as the structure of the donor star and its re- 
sponse to mass loss can be described using relatively simple 
prescriptions; however, this is not possible in general. The clear 
and widespread realization that intermediate-mass donor stars 
can stably transfer matter to a neutron star (NS) accretor came 
largely as a result of recent calculations that utilized full stellar 
evolution codes (e.g., Tauris & Savonije 1999; Podsiadlowski 
& Rappaport 2000; Tauris et al. 2000; Kolb et al. 2000; Podsi- 
adlowski et al. 2002). Only with such codes can the evolution 
of the donor be followed realistically during the rapid phase 



of thermal timescale mass transfer that characterizes the early 
evolution of intermediate-mass X-ray binaries (IMXBs). 

Podsiadlowski et al. (2002, hereafter, Paper I) is devoted to 
a systematic evolutionary study of L/IMXBs, wherein we de- 
scribe a library of 100 evolutionary sequences computed with 
a standard Henyey-type stellar structure code. This library has 
now been expanded to 144 sequences, covering initial orbital 
periods from 2 hours to 100 days and initial donor masses from 
0.3 to 1 Mq. The library is intended to provide a fairly complete 
mapping of the initial conditions that are likely to be encoun- 
tered in a population synthesis study of L/IMXBs. 

Here we extend the work in Paper I by combining our library 
with a detailed Monte Carlo binary population synthesis (BPS) 
code for L/IMXBs. The code includes standard assumptions 
for the population of massive primordial binaries, reasonable 
analytic prescriptions to describe both stable and dynamically 
unstable mass transfer prior to the supernova (SN) explosion, 
and NS kicks. Similar codes are described in Portegies Zwart 
& Verbunt (1996) and Belczynski et al. (2002). 

We undertake a limited exploration of the set of free parame- 
ters that enter the BPS calculation used to generate the incipient 
X-ray binaries. Probably the most important parameters in the 
BPS study are the mean NS kick speed and the envelope bind- 
ing energy that enters the prescription for common-envelope 
evolution. For reasonable variation of these two quantities, the 
formation probability of L/IMXBs ranges over two orders of 
magnitude. 

For each incipient L/IMXB that emerges from the BPS cal- 
culation, we find an initial model in our library with the closest 
matching orbital period and donor mass. For the ensemble of 
selected sequences, we apply a temporal weighting scheme to 
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calculate the distributions of potentially observable quantities 
at the current epoch. This is the first paper where such distribu- 
tions have been computed for L/IMXBs, and it is now possible 
to directly compare population models and the statistics of ob- 
served systems. 

The paper is organized as follows. In § 2, we briefly describe 
our BPS code, highlighting the important uncertainties and the 
associated free parameters. The population of incipient X-ray 
binaries that emerges from the BPS calculation is discussed in 
§ 3. Key results of this study are presented in § 4, where we 
show distributions of various quantities at the current epoch and 
make rough comparisons with the observational data. Finally, 
in § 5, the most important results of our investigation are listed, 
along with suggestions for how this work may be extended and 
improved. 

2. MASSIVE BINARY POPULATION SYNTHESIS 

The formation of a NS in a binary system involves three main 
evolutionary steps: (1) the formation of a primordial binary, 
where the initially more massive component (the primary) has 
a mass > 8M Q , (2) a phase of mass transfer from the primary 
to the secondary (the initially less massive component), and 
(3) the subsequent SN explosion of the primary's hydrogen- 
exhausted core and the formation of the NS. Our Monte Carlo 
BPS code utilizes a set of analytic prescriptions to describe each 
of these steps. A brief overview of the important elements of 
our code is given below; an expanded account is provided in 
Pfahl et al. (2002). 

2.1. Primordial Binaries 

We construct each primordial binary by selecting the compo- 
nent masses and orbital parameters from the following distribu- 
tion functions. 

Primary Mass. — The initial primary mass, Mu, is chosen from 
a power-law initial mass function, p(Mu) oc M^f. We use a fixed 
value of x = 2.5 for massive stars (e.g., Miller & Scalo 1979; 
Scalo 1986; Kroupa et al. 1993). Primary masses are restricted 
to the range Mu = 8-25 M Q , and we assume that the primary 
is always the NS progenitor (see, however, Podsiadlowski et al. 
1992; Pols 1994; Wellstein et al. 2001). 

Secondary Mass. — The initial secondary mass, Mu, is cho- 
sen from a distribution in mass ratios, p(qf) oc qj, where qi = 
M2i/M\i < 1 . Strongly motivated by the work of Garmany et al. 
(1980), we prefer a flat distribution (y = 0), but we also consider 
y = -1 and y = 1. 

Eccentricity. — Without much loss in generality, we take the pri- 
mordial binary orbits to be circular. This assumption is dis- 
cussed in Pfahl et al. (2002). 

Semimajor Axis. — The initial orbital separation, a,-, is drawn 
from a distribution that is uniform in loga,- (e.g., Abt & Levy 
1978). We determine the minimum value of a,- for each system 
by demanding that neither star overflows its Roche lobe on the 
main sequence. The upper limit is somewhat arbitrary, but here 
is taken to be 10 3 AU. 

2.2. Mass Transfer 

If a, < 5-10 AU, the primary will grow to fill its Roche lobe 
at some point during its evolution. The subsequent phase of 
mass transfer is of crucial importance in determining what types 



of NS binaries are ultimately produced. It is common prac- 
tice to distinguish among three main evolutionary phases of the 
primary at the onset of mass transfer (Kippenhahn & Weigert 
1967; Podsiadlowski et al. 1992). Case A evolution corre- 
sponds to core hydrogen burning, case B refers to the shell 
hydrogen-burning phase, but prior to central helium ignition, 
and case C evolution begins after core helium burning. It is 
quite improbable for mass transfer to begin during core helium 
burning, and we thus neglect this possibility (see Pfahl et al. 
2002). We refer to as case D the large fraction of wide binaries 
that remain detached prior to the SN explosion of the primary. 
Using the distributions and standard-model parameters given 
above, as well as the treatment of stellar winds discussed be- 
low, we find that cases A, B, C, and D comprise roughly 5%, 
30%, 15%, and 50%, respectively, of the primordial binary pop- 
ulation. In order to determine which case each binary falls into, 
we use the single-star evolution fitting formulae of Hurley et al. 
(2000). 

Mass transfer from the primary to the secondary may be sta- 
ble or dynamically unstable, depending mainly on the binary 
mass ratio and evolutionary state of the primary when it fills its 
Roche lobe. In our population study of L/IMXBs, we consider 
only cases B and C mass transfer. Case A mass transfer ac- 
counts for only a small fraction of the binaries and, furthermore, 
most likely leads to the merger of the two stars following a con- 
tact phase (e.g., Wellstein et al. 2001). Most case D systems are 
disrupted due to the SN explosion of the primary if NS kicks 
are significant. We do not consider case D systems that survive 
the SN; for a discussion of the products that may emerge from 
this evolutionary channel, see Kalogera (1998) and Willems & 
Kolb (2002). 

Cases B and C are divided in early (B e or C e ) and late (B; or 
Q) phases, if the primary has an envelope that is mostly radia- 
tive or deeply convective, respectively. We assume that mass 
transfer is stable, though non-conservative, for cases B e and C e 
if the mass ratio, after any wind mass loss has occurred, is >q c , 
where q c is some critical mass ratio. We adopt a fixed value of 
q c = 0.5 in our study (e.g., Wellstein et al. 2001). If the primary 
has a deep convective envelope when it fills its Roche lobe, 
mass transfer is dynamically unstable and a common-envelope 
(CE) phase ensues, which results in either a very compact bi- 
nary or a merger. 

A single star of mass >15M Q may lose >30% of its mass 
in a stellar wind on the asymptotic giant branch (AGB). For 
stars of mass <25M , only <5% of the mass is lost on the 
main sequence. We suppose that the wind from the primary in 
a binary system takes with it the specific orbital angular mo- 
mentum of the star. In response to the AGB winds, the Roche 
lobe of the primary expands and may overtake the expansion of 
the star, making Roche-lobe overflow and case Q mass transfer 
impossible. We have included the effects of stellar winds only 
for initial primary masses >13M Q , on both the main sequence 
and the AGB; our procedure is similar to the one adopted by 
Podsiadlowski et al. (2003). For this range of masses, core he- 
lium burning begins while the star is in the Hertzsprung gap, 
and there is no decrease in the stellar radius. Primaries of 
mass Mu < 13M Q experience moderate wind mass loss dur- 
ing core helium burning following evolution through the first 
giant branch, but the stellar radius decreases after helium ig- 
nition, precluding Roche-lobe overflow during this phase. We 
note that the mass that separates the two behaviors just men- 
tioned is actually quite uncertain (e.g., Langer & Maeder 1995), 
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and may be as large as ^20 M . 

If a merger is avoided, it is reasonable to suppose that the 
primary loses its entire envelope, leaving only its hydrogen- 
exhausted core, irrespective of whether mass transfer is stable 
or dynamically unstable. Following case B mass transfer, the 
mass of the helium core is given approximately by (Hurley et al. 
2000) 

M d ~ 0.1M, 1 , 35 . (1) 

We neglect the relatively small amount of wind mass loss on the 
main sequence and use the initial mass of the primary to com- 
pute the initial core mass M c ,- before it too loses mass in a wind. 
The mass of the exposed core may be larger by ^0.5-1 M Q af- 
ter case C mass transfer, as a result of shell nuclear burning. 
A helium core will ultimately yield a NS remnant if its mass 
is >2Mq (Habets 1986b). Equation (1) gives a primary mass 
threshold for NS formation of ^9M Q 

For our chosen maximum primary mass of 25 M Q , the corre- 
sponding core mass is ^8M Q . A nascent helium star of mass 
3-8 Mq that is exposed following case B mass transfer may lose 
10-30% of its mass in a wind before the SN (e.g., Brown et al. 
2001; Pols & Dewi 2002). The final core mass is related to the 
initial helium star mass by the approximate formula (see Fig. 1 
of Pols & Dewi 2002) 

M cf ~lAM 2 J\ (2) 

for M ci > 3M Q . Following case C mass transfer, the core of 
the primary has already undergone helium burning, and there is 
insufficient time for winds to significantly reduce its mass. 

If M C i ^ 3M following case B mass transfer, we may 
safely neglect winds, but such helium stars may expand to gi- 
ant dimensions following core helium burning, often initiating 
a phase of so-called case BB mass transfer to the secondary 
(De Greve & De Loore 1977; Delgado & Thomas 1981; Ha- 
bets 1986a). We do not attempt to model this evolution in de- 
tail, but simply assume that O.5M is transferred conservatively 
from the primary's core to the secondary. In case BB systems 
where M 2 i < M„, such as when the secondary is of low- or 
intermediate-mass, the mass transfer may proceed on the ther- 
mal timescale of the core, and the evolution may be quite com- 
plicated (e.g., Dewi et al. 2002; Ivanova et al. 2003). However, 
any reasonable treatment of case BB mass transfer is not likely 
to change our results for L/IMXBs substantially, 

Stable mass transfer from the primary to the secondary is 
treated analytically as follows. We assume that the secondary 
accretes a fraction (3 of the material lost from the primary dur- 
ing Roche-lobe overflow. The complementary mass fraction, 
1-/3, escapes the system with specific angular momentum a, 
in units of the orbital angular momentum per unit reduced mass. 
We use constant values of a = 1.5, characteristic of mass loss 
through the L2 point, and (3 = 0.75. The final orbital separa- 
tion is then given by the generic equation (Podsiadlowski et al. 
1992) 
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where 



Ci =2a(l-/?)-2 

C 2 =-2a(l-/?)//3-2. (4) 

Here the subscript 'RLO' denotes stable Roche-lobe overflow, 
and unprimed and primed quantities denote parameters at the 
onset and termination of mass transfer, respectively. 



It is easily verified that, in our simulations, the minimum 
secondary mass resulting from stable mass transfer is roughly 
q c SM & + /36Mq, where 8M Q is the minimum primary mass 
and 6M Q is the corresponding envelope mass shed during mass 
transfer. For our chosen values of q c = 0.5 and j3 = 0.75, this 
minimum mass is 8.5M Q , considerably larger than the maxi- 
mum initial donor mass of ^4M for which an IMXB under- 
goes stable mass transfer (see Paper I). Thus, in our work, all 
incipient L/IMXBs are the products of dynamically unstable 
mass transfer. 

We use the conventional energy relation to describe the dy- 
namical spiral-in during a CE phase (e.g., Webbink 1984). The 
ratio of the final to the initial orbital separation is given by the 
generic equation 

-l 



M e M 2 



M 2 + ■ 



2M„ 



(5) 



CE M i \ AcE^CEru, 
where ru is the Roche-lobe radius of the primary in units of the 
orbital separation, M e and M c are the masses of the primary's 
envelope and core, respectively, Ace parameterizes the struc- 
ture of the envelope, and tjce is the fraction of orbital binding 
energy that goes into dissipating the envelope. Although Ace 
and ?7ce appear only in the product Ace Vce, we suppose that 
tjce = 1 and fix Ace in the range 0.1-0.5 for each simulation. 
In reality, Ace changes as a star evolves, typically decreasing 
to <0. 1-0.2 as the massive star passes through the Hertzsprung 
gap, and increasing to ^0.1-0.4 as the star ascends the AGB 
(Dewi & Tauris 2000, 2001). Generally, Ace decreases for all 
stellar radii as the mass of the star is increased. Equation (5) 
gives typical shrinkage factors of (a'/a)cE ~ 0.01 for systems 
with low- and intermediate-mass secondaries. 

A sufficient condition for the merger of the primary and sec- 
ondary is that the main-sequence secondary overfills its Roche 
lobe for the calculated post-CE orbital separation. Therefore, 
the minimum separation for surviving systems must be larger 
than several solar radii, corresponding to initial orbital separa- 
tions greater than several hundred solar radii. It turns out that 
the majority of the dynamically unstable case B e and C e sys- 
tems merge following the CE, and that in most systems that 
survive the CE, the primary is a convective red supergiant (case 
B; or Q) at the onset of mass transfer. 

2.3. Supernova Explosion 

After the exposed core of the primary consumes its remain- 
ing nuclear fuel, it explodes as a Type lb or Ic SN and leaves a 
NS remnant. We take the initial NS mass to be 1.4M Q . Impul- 
sive mass loss and the NS kick strongly perturb the binary and 
may cause its disruption. Mass loss is especially significant, 
since the mass ejected may be comparable to or greater than 
the secondary mass. Some important insights can be obtained 
rather simply by neglecting NS kicks. 

It is straightforward to show (e.g., Blaauw 1961; Boersma 
1961) that for a circular pre-SN orbit and vanishing kicks the 
eccentricity, esN, after the SN is simply 

M C -M NS 

esN = — — 7T— , (6) 
M 2 +M NS 

where M N s is the mass of the NS, and M c is the pre-explosion 
core mass. The system is unbound when M c > M 2 + 2M^s,- Let- 
ting M 2 = 1M , we see that for M c > 3.8M Q (Mi > 15M ) dis- 
ruption of the binary is guaranteed. A kick of appropriate mag- 
nitude and direction is then required in order to keep the sys- 
tem bound. For intermediate-mass secondaries, a wider range 
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of pre-SN core masses is permitted when kicks are neglected. 
Equation (6) gives values of esN ^ 0.5 for some typical core and 
secondary masses. Following the CE phase, the core of the pri- 
mary will have orbital speeds about the binary center-of-mass 
(CM) of v c > 200km s" 1 . After the SN, mass loss alone gives 
the binary CM a speed of vcm = <?snv c , so tnat eccentricities of 
esN ^ 0.5 correspond to large post-SN systemic speeds of order 
lOOkms" 1 . 

When NS kicks are considered in addition to SN mass, a 
larger fraction of systems are disrupted, and those binaries that 
do remain bound will have larger CM speeds. We utilize a 
Maxwellian distribution in kick speeds, 

p(v k )=J-^exp(-v 2 k /2a 2 k ), (7) 

where the directions of the kicks are distributed isotropically. 
Dispersions of cr* ~ 100-200 km s" 1 are reasonably consistent 
with the data on pulsar proper motions (e.g., Hansen & Phin- 
ney 1997; Arzoumanian et al. 2002). However, neither the 
functional form of the kick distribution nor the mean are very 
well constrained. In our study, we consider = 50, 100, and 
200km s -1 . The post-SN orbital parameters are calculated using 
the formalism described in Appendix B of Pfahl et al. (2002). 

Significant SN mass loss and large NS kicks yield bound 
post-SN binaries with high eccentricities. Given the post-SN 
eccentricity, we check if the radius of the secondary is >10% 
larger than its tidal radius at periastron. If this occurs, we as- 
sume that the NS immediately spirals into the envelope of the 
unevolved secondary and do not consider the system further. 
Our 10% overflow restriction allows for the possibility that tidal 
circularization and perhaps some mass loss will prevent the ob- 
jects from merging. The details of eccentric binary evolution 
with mass transfer is well beyond the scope of this investiga- 
tion. However, we note in passing that Cir X-l is most likely 
a young, possibly intermediate-mass, X-ray binary undergoing 
episodic mass transfer as a result of having a highly eccentric 
(e SN > 0.8) orbit (Shirey 1998). 

If the coalescence of the NS and secondary is avoided, we 
neglect mass loss from the system and assume, rather simplisti- 
cally, that the binary circularizes while conserving orbital angu- 
lar momentum. The final orbital separation is then «sn(1 -<?g N ), 
where «sn is the semimajor axis after the SN. Moreover, we 
assume that the secondary rotates synchronously with the cir- 
cularized orbit. 

3. INCIPIENT X-RAY BINARIES 

The output of the population synthesis calculation is a set of 
circular binaries, each identified by their orbital period and the 
mass of the secondary. In order to select initial models from 
our library of L/IMXB evolutionary sequences, we require the 
orbital period at which the secondary first fills its Roche lobe. 
Hereafter, we denote the donor and accretor (NS) masses in 
units of Mq by Md and M a , respectively. 

Orbital angular momentum losses via gravitational radiation 
(GR) and magnetic braking (MB) may cause the binary sep- 
aration to decrease substantially prior to mass transfer. The 
timescale for orbital shrinkage due to GR is 

/ p \ 8/3 

TGR-llOGyr^a+^/M ^g-J , (8) 

where q = Md/M a , and M a =l .4. We utilize the same standard 
MB formula as in Paper I (Verbunt & Zwaan 1981; Rappaport 



et al. 1983), with a characteristic timescale given by 

/ p . \ 10/3 

tub 3.5Gyr ^(1 +qT 1/3 R d 4 i J , (9) 

where Rd is the radius of the secondary in solar units, and we 
introduce 7/mb to parameterize the strength of MB prior to mass 
transfer. Only when the secondary is tidally locked to the orbit 
does MB extract orbital angular momentum. In an ad hoc way, 
we let i]mb = or 1 to model the cases where the tidal coupling 
prior to Roche-lobe overflow is very weak or very strong, re- 
spectively. We apply MB only if 0.4 < Md < 1.5; less massive 
stars are fully convective and may not undergo MB (Rappaport 
et al. 1983), while more massive stars have radiative envelopes 
and may not undergo MB. When the main sequence lifetime of 
a low-mass star exceeds the ~13 Gyr age of the Galaxy, MB or 
GR are required to drive the system into contact. Intermediate- 
mass secondaries have nuclear lifetimes of ~0.1-3 Gyr, and so 
do not suffer from this obstacle. 

The single-star evolution code of Hurley et al. (2000) is used 
to follow the radial evolution of the secondary as MB and GR 
shrink the orbit. The metallicity is set to the solar value of 
Z = 0.02. We neglect the evolutionary time prior to the SN of 
the primary and assume that the secondary is on the ZAMS at 
the time the NS is formed. The age of the secondary when it 
fills its Roche lobe is here referred to as the "lag time," denoted 
by fi a g. The lag time thus approximates the time between the 
formation of the primordial binary and the onset of the X-ray 
binary phase. We take the age of the Galaxy to be 13 Gyr, and 
only accept as incipient L/IMXBs those for which fi ag < 1 3 Gyr. 

Figure 1 shows an illustrative set of distributions of binary 
parameters and systemic speeds for the incipient L/IMXBs. 
Overlayed on the scatter plots are the initial models (open cir- 
cles) from our library of L/IMXBs evolutionary sequences. We 
applied fixed values of y = 0.0, at = 200km s" 1 , i] M B = L and 
considered CE structure parameters of Ace = {0.1,0.3,0.5}. 
For each parameter set shown, the distribution of initial sec- 
ondary masses increases with mass. This distribution is mod- 
ified if we adopt a different distribution of mass ratios for the 
components of the primordial binaries. Nonetheless, the gen- 
eral statistical importance of initially intermediate-mass sec- 
ondaries is clear. The high concentration of systems with 
P orb < 0.5 d and M d < 1 .5 when Ace = 0.3 and 0.5 is due to MB 
and GR. The absence of LMXBs with evolved secondaries of 
mass Md < 1 for f or b > 0.5 d when Ace = 0.3 and 0.5 results 
from the demand that fi ag < 13 Gyr. More low-mass donors 
survive this cut if the metallicity is reduced, since the main- 
sequence lifetime decreases with decreasing metallicity. Note 
that for Ace = 0.1 no systems with low-mass donors survive the 
CE, and that the maximum period of the incipient IMXBs is 
only ^3 day. When Ace = 0.1, the small fraction of incipient 
IMXBs with 2.5 < Md < 4 is not well sampled by our grid of 
evolutionary sequences, but we expect that a denser sampling 
will not greatly change our results. 

It is somewhat interesting that systems with fully convective 
secondaries of mass Md < 0.4 can be driven into contact by GR 
within the lifetime of the Galaxy (eq. [8]). A handful of such bi- 
naries are present in Fig. 1 for Ace = 0.3 and 0.5. Mass transfer 
will be driven by GR at rates of <10" 10 M Q yr" 1 , from initial pe- 
riods of ^2-3 hr to a minimum period of ~ 1 .5 hr, at which point 
Md ~ 0.05. The secondary and orbit then expand in response 
to further mass transfer. For a physical discussion of LMXB 
evolution with very low-mass donors and the period minimum, 
see Nelson et al. (1986). The companion mass (M d > 0.04) 
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and orbital period (P or b ~ 2hr) of the first-discovered millisec- 
ond X-ray pulsar SAX J1808.4-3658 (Wijnands & van der Klis 
1998; Chakrabarty & Morgan 1998) are consistent with it hav- 
ing formed and evolved in the way just described. However, 
given the uncertainty in the binary inclination, a model wherein 
the initial donor mass takes a more typical value of ~1M Q is 
also consistent with the present system parameters (L. Nelson 
& S. Rappaport 2003, in preparation). 

We define the formation efficiency as T(<Md) = M nc (< 
Md) /Npb, where Ni nc (<Md) is the synthesized number of incip- 
ient L/IMXBs with donor masses less than Md, and Npb is the 
number of primordial binaries used in the simulation. The ma- 
jority of systems with Md > 4 undergo dynamical mass transfer 
shortly after Roche-lobe overflow (see Paper I). Therefore, we 
take !F(<4) to be the total formation efficiency of LMXBs and 
IMXBs. We somewhat arbitrarily define incipient LMXBs as 
systems with initial donor masses of Md < 1.5. The present 
rate of core-collapse SNe in the Milky Way is ~ 10~ 2 yr _1 
(e.g., Cappellaro et al. 1999). Here we adopt the simplifying 
assumption that the average star formation rate has been con- 
stant for 13 Gyr and take 7?-sn to be the approximate formation 
rate of massive primordial binaries. It follows that the current 
Galactic birthrate of L/IMXBs is B(<M d ) = H<M d )TZsN- 

Table 1 lists formation efficiencies for a modest number of 
different parameter sets in the BPS study, where we vary only 
y, Ace, &k, and ?7mb- We used 10 7 primordial binaries in each 
simulation. For the parameter sets explored, log J r (<4) varies 
from roughly -4.1 to —1.7, with corresponding birth rates of 
B«4) - 10- 6 -10- 4 yr _1 . Incipient LMXBs are formed with 
systematically lower efficiencies of log.F(<1.5) < -2.4, and 
birthrates of S(<1.5) < 4 x 10 _5 yr _1 - The ran g e of LMXB 
birthrates is reasonably consistent with the results of, e.g., 
Portegies Zwart & Verbunt (1996) and Kalogera & Webbink 
(1998), each of whom used somewhat different assumptions 
than the ones adopted here. The formation efficiencies for 
Md < 0.4 are also listed in Table 1, with typical values of 
logJF(<0.4) < -4. If SAX J1808.4-3658 evolved in the way 
described above, then Ace ~ 0.5 if favored, since this value 
yields the largest birthrates of B(< 1 .5) ~ 10" 6 yr" 1 . 

4. POPULATION AT THE CURRENT EPOCH 

In order to compare our theoretical population synthesis re- 
sults and the statistics of observed systems, it is necessary to 
follow the evolution of each L/IMXB that we generate. For 
each synthesized system, a close match is found in our library 
of evolutionary sequences. Distributions at the current epoch of 
observable quantities are computed by appropriately weighting 
each evolutionary sequence. We now elaborate on these points 
and discuss several key results, as well as make very rudimen- 
tary comparisons between our results and the observations. 

4.1. Selection from the Library 

Each incipient L/IMXB is characterized by the donor mass 
and the circularized orbital period at which the star first fills its 
Roche lobe. For any given initial P or b and Md, we select a model 
from our L/IMXB evolutionary library by first identifying the 
subset of sequences in the library with the closest initial donor 
mass, and then finding the one sequence in this subset with the 
closest initial orbital period. Incipient LMXBs with Md < 0.3 
are evolved by selecting the one sequence with an initial donor 
mass of 0.3M Q and considering only the part of the sequence 
for which the donor mass is <Md- This is reasonable, since 



such stars will initially follow the approximate mass-radius re- 
lation, Rd oc M ( ^ & , for low-mass stars in thermal equilibrium. 
If the library sequences are labeled with some integer index i, 
we may define the formation efficiency J 7 , of sequence i as the 
number of times that sequence is selected, divided by Npb- 

4.2. Weighting Procedure 

For each selected library model, the entire evolutionary se- 
quence contributes to the distributions at the current epoch of 
^orb, Mj, andM fl , the accretion rate onto the NS. This is done as 
follows. Let Q be the quantity of interest. The evolutionary data 
file for each sequence gives Q as a function of the time ?mt since 
the onset of mass transfer to the NS. For some small interval of 
time 5t M T, Q varies over a small range (£>', Q! + 6Q). The prob- 
ability that there is an identical system that is presently in this 
particular state somewhere in the Galaxy is -F/T^sn^pb^mt, 
where f PB = — (?mt + flag) is the formation time of the primor- 
dial binary, and now is taken to be at time t = 0. The total 
number of L/IMXBs present in the Galaxy with Q in the bin 
(2o,2o + A0 is ^-T^sn^peO^mt- Here the symbolic sum 
is over all sequences i and all evolutionary times for which Q 
lies in the chosen bin. While the above procedure is valid for 
a variable star formation rate, we limit our study to a constant 
massive binary formation rate of TZsn = 10" 2 yr _1 (see § 3). 

4.3. Current-Epoch Distributions 

Figure 2 shows theoretical distributions (solid histograms) at 
the current epoch of Md, P ox b, and M a , where we have used 
illustrative values of y = 0.0, a k = 200kms _1 , and rj MB = 1. 
The left and right panels are for Ace = 0.1 and 0.5, respec- 
tively. Hatched regions indicate systems that are persistent X- 
ray sources (+45°) or transient (-45°) (see below). The solid 
histogram that encloses the hatched regions shows the total 
number of systems. The dotted and thick, solid histograms in 
the right panel are explained in § 4.4 and § 5.2. 

We decide if a system is persistent or transient according to 
the standard disk instability model (e.g., Cannizzo et al. 1982), 
wherein the accretion is transient if the X-ray irradiation tem- 
perature at the disk edge is T m > 10 4 K (see van Paradijs 1996, 
and references therein), a characteristic hydrogen ionization 
temperature. The irradiation temperature at a radius R in the 
disk is given by (e.g., de Jong et al. 1996) 

r -=4^f (1 -^- 1) ' (10) 

where osb is the Stefan-Boltzmann constant, H is the disk scale 
height, 7 is the X-ray albedo of the disk, and £ = dlnH/dlnR. 
Following, e.g., King et al. (1996), we take the outer disk ra- 
dius to be 70% of the Roche-lobe radius of the NS, and adopt 
values of £ = 9/7 (Vrtilek et al. 1990), 7 = 0.9, and H/R = 0.2 
at the disk edge (de Jong et al. 1996). Fixing the NS mass and 
radius at 1.4M Q and 10 km, respectively, we find the follow- 
ing expression for the critical mass-transfer rate below which 
an L/IMXB is transient: 

lO-i°M vr-i ~ (1+g) MihrJ ' (U) 

where is the Roche-lobe radius for the NS in units of the 
orbital separation. The value of M^crit is not certain to within 
a factor of at least a few, and so our numbers and distributions 
for transient and persistent systems should be considered cau- 
tiously. 
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Table 1 

Selected results from our population study 



Model 
Number 




A CE b 


»?MB C 




logJF(<4) e 


logJFKl.5) 


log T{<0A) 


logJV x f 


/sp s 


<TX> 


1 


0.0 


0.5 


1.0 


200 


-2.17 


-2.89 


-3.99 


4.52 


0.36 


5.03 


2 


0.0 


0.5 


1.0 


100 


-1.86 


-2.60 


-3.77 


4.76 


0.35 


4.28 


3 


0.0 


0.5 


1.0 


50 


-1.70 


-2.49 


-3.84 


4.83 


0.35 


3.53 


4 


0.0 


0.3 


1.0 


200 


-2.38 


-3.22 


-5.09 


4.37 


0.31 


5.73 


5 


0.0 


0.3 


1.0 


100 


-2.06 


-2.91 


-5.00 


4.62 


0.30 


4.91 


6 


0.0 


0.3 


1.0 


50 


-1.90 


-2.74 


-5.09 


4.66 


0.29 


3.77 


7 


0.0 


0.1 


1.0 


200 


-3.87 


< -7 


< -7 


2.84 


0.30 


5.23 


8 


0.0 


0.1 


1.0 


100 


-3.56 


< -7 


< -7 


3.14 


0.29 


5.12 


9 


0.0 


0.1 


1.0 


50 


-3.41 


< -7 


< -7 


3.30 


0.23 


5.33 


10 .... 


0.0 


0.5 


0.0 


200 


-2.20 


-3.04 


-3.99 


4.49 


0.31 


4.92 


11 .... 


0.0 


0.3 


0.0 


200 


-2.40 


-3.33 


-5.09 


4.35 


0.27 


5.73 


12 .... 


0.0 


0.1 


0.0 


200 


-3.87 


< -7 


< -7 


2.84 


0.30 


5.22 


13 .... 


-1.0 


0.5 


1.0 


200 


-2.02 


-2.38 


-3.03 


4.63 


0.46 


4.63 


14 .... 


-1.0 


0.3 


1.0 


200 


-2.32 


-2.81 


-4.30 


4.43 


0.38 


5.75 


15 .... 


-1.0 


0.1 


1.0 


200 


-4.06 


< -7 


< -7 


2.66 


0.32 


5.36 


16 .... 


1.0 


0.5 


1.0 


200 


-2.56 


-3.71 


-5.52 


4.06 


0.29 


4.30 


17 .... 


1.0 


0.3 


1.0 


200 


-2.71 


-3.98 


-6.00 


4.00 


0.26 


5.25 


18 .... 


1.0 


0.1 


1.0 


200 


-4.06 


< -7 


< -7 


2.62 


0.29 


5.14 



a Exponent for the mass-ratio distribution. 

b Binding-energy parameter for the common envelope. 

c Strength parameter for magnetic braking. 

d Dispersion (in kms -1 ) of the Maxwellian kick distribution. 

e To obtain the approximate birthrates, multiply T by 7£sn = 10~ 2 yr -1 . 

f Total number of systems at the current epoch with M a = 10 -10 -10 -8 M & yr -1 . 

g Fraction of short-period (P or b < 2hr) systems at the current epoch. 

h Mean X-ray lifetime in units of 10 8 yr for M a = 10~ 10 -10~ 8 M & yr _1 . 



In Fig. 2, the shapes and extents of the theoretical distribu- 
tions for Ace = 0.1 and 0.5 are quite similar. The most notable 
difference between the two cases is the total number of sys- 
tems: ~T0 3 when Ace = 0.1 and ~10 5 when Ace = 0.5. Also 
prominent in both panels is the very large number of systems 
with P rb 2hr, the majority of which are transient according 
to eq. (11). Note finally that, while IMXBs are much more 
favorably produced than LMXBs, very few systems at the cur- 
rent epoch have donors of mass > 1 M@. This is simply because 
the initial phase of thermal-timescale mass transfer in IMXBs, 
where a large fraction of the secondary mass is removed, is rel- 
atively short-lived (see, e.g., Fig. 2a of Paper I). 

4.4. Comparisons with Observation 

Small number statistics, observational selection effects, and 
sample incompleteness are all serious issues for the observed 
LMXBs. Out of ^140 Galactic X-ray sources classified as 
LMXBs with NS accretors, there are only ^40 systems with 
measured orbital periods and a handful with estimated sec- 
ondary masses. Usually, the nondetection of the donor at opti- 
cal wavelengths is taken to mean that it is of low mass. The esti- 
mation of X-ray luminosities among the observed LMXBs that 
lie outside of globular clusters is complicated by very poor dis- 
tance estimates. Furthermore, for X-ray luminosities less than 
L x ~ 10 36 ergss _1 (M a ~ 10" 10 M Q yr" 1 ), the observed sample 
may be quite incomplete. Even in light of these problems, some 
rough comparisons between our theoretical results and the ob- 
servations are illuminating. The most apparent quantitative dis- 
crepancy is in the total numbers of systems. 



We define N x to be the total number of L/IMXBs at the cur- 
rent epoch — both persistent and transient — with secular accre- 
tion rates of M a > 10~ 10 M Q yr -1 . Thus, iVx quantifies the num- 
ber of luminous X-ray binaries that would be observable over 
a large fraction of the Galactic volume. Table 1 lists values of 
logA^x for different parameter sets, from which we see that Nx 
ranges from ^3 to >400 times the the total of ^140 systems 
observed in the Galaxy. We elaborate in the next section on this 
rather severe overproduction problem. 

In Fig. 2, for Ace = 0.5, we have overlayed (dotted histogram) 
on the middle panel the distribution of 37 measured orbital pe- 
riods for LMXBs outside of globular clusters Liu et al. (2001), 
which has been multiplied by a factor of 3000 for compari- 
son purposes. With such small numbers (1-9) per bin, mean- 
ingful comparisons with our theoretical distribution are diffi- 
cult. It may be that we are underproducing systems with P or b > 
lOd and overproducing short-period LMXBs with P or b < 0.1 d. 
Our simulations indicate that most of the short-period systems 
are transient, which may aid in explaining a possible discrep- 
ancy. Furthermore, binaries with P orD < 2 hr are driven by GR 
with characteristically low X-ray luminosities of Lx < 10 35 - 
10 36 ergss _1 , making them more difficult to discover. We com- 
pute the fraction, fsp, of short-period (P orD < 2hr) LMXBs at 
the current epoch, and find that typically fsp ~ 0.3 (Table 1), 
which in itself is not in serious conflict with observations, for 
reasons already mentioned. 

In our simulations, we actually underestimate the number of 
short-period (P orD < 2hr) LMXBs at the current epoch. These 
binaries evolve from systems with initial orbital periods below 



Low- and Intermediate-Mass X-ray Binaries 



7 



o ■ " oo d" o" 6' <>' 

o »oo ° ° « 9 ' 
o, ° > a o 

o o-4^o o. ■<>' O -C. 6 " 



0.5 

200 km s- 



1 2 3 4 5 6 
Secondary Mass (M ) 




-1 -0.5 0.5 1 1.5 
log Period (days) 



1 2 3 4 5 6 7 
Secondary Mass (M e ) 





50 100 150 200 250 300 

Systemic Speed (km s _1 ) 




1 2 3 4 5 6 
Secondary Mass (M ) 




1 2 3 4 5 6 
Secondary Mass (M Q ) 




<u 
OS 



-1 -0.5 0.5 1 1.5 
log Period (days) 




50 100 150 200 250 300 
Systemic Speed (km s _1 ) 



o ° o ° 
o ° 

o °° °°° 

o°°° o 

o° °° o o o 

"o. ° o o o o 

°oo °o o o o o 
ooo q00 o o o o 

og° 
8 



0.1 

200 km s" 1 

J I I I l_ 



1 2 3 4 5 6 
Secondary Mass (M ) 




1 2 3 4 5 
Secondary Mass (M s ) 




<u 
06 



-0.5 0.5 1 
log Period (days) 




50 100 150 200 250 300 
Systemic Speed (km s _1 ) 



FIG. 1. — Scatter plots and histograms for incipient L/IMXBs. For these simulations we used 500,000 primordial binaries and parameter values of y = 0.0, 
o> = 200kms~' , ?/MB = 1. and Ace = {0.1,0.3,0.5}. The clustering of systems with P or b < 0.5 d and Mj < 1.5 for Ace = 0.3 and 0.5 is due to the effects of magnetic 
braking and gravitational radiation. The histograms also show the distributions of systemic speeds that result from the SNe. For reference, we have overlayed on the 
scatter plots the initial models in our library of evolutionary sequences (open circles). 
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FIG. 2. — Distributions at the current epoch of the orbital period, donor mass, and mass accretion rate onto the NS, for parameters y = 0.0, at = 200 kms -1 , 
r?MB = 1. and Ace =0.1 (left) and 0.5 (right). We have assumed a constant massive binary formation rate of 10~ 2 yr~' . Note the difference in scale between the left 
and right figures. The hatched regions indicate persistent (+45°) and transient (-45°) X-ray sources, and the enclosing solid histogram gives the sum of these two 
populations. Overlayed (dotted histogram) on the theoretical period distribution in the figure on the right is the rescaled distribution of 37 measured periods (Liu 
et al. 2001) among ~140 observed LMXBs in the Galactic plane. The thick, solid distribution overlayed on the bottom panel of the right figure illustrates how the 
inclusion of X-ray irradiation effects might change our theoretical M a distribution. 



the bifurcation period of ~18 hr (see Paper I and references 
therein), and reach minimum periods of ^lOmin to ~1.5 hr. 
However, for technical reasons, our X-ray binary calculations 
are terminated not long after the minimum period is reached. 
The subsequent evolution is driven by GR at low rates, but may 
last for billions of years. This makes the discrepancy with the 
observations somewhat worse. On the other hand, we note that 
3 ultracompact binary systems have been discovered within the 
past year alone, so the discovery probability for these systems 
may be increasing. 

Extrapolating the data available for 16 Galactic LMXBs 
monitored by RXTEI 'ASM, Grimm et al. (2002) have attempted 
to compute a cumulative X-ray luminosity distribution, cor- 
rected for the fraction of the Galactic volume observable by 
the ASM. For luminosities of >10 36 ergs s" 1 , we expect that the 
observed sample should be reasonably complete. As a point of 
comparison, we have chosen to compute the ratio of the number 
of LMXBs with X-ray luminosities of L x = 10 36 -10 37 ergss _1 
to the number with >10 37 ergss _1 . The results of Grimm et al. 
(2002) indicate that ratio is ~1.5. Our theoretical cumulative 
distribution of mass accretion rates is shown as the thin line in 
Fig. 4; the thick line will be discussed in § 5.2. For the two 
luminosity ranges given above, we find a larger number ratio of 
~5.3. It is not clear if this particular discrepancy between our 
theoretical results and the observations is especially significant. 

5. DISCUSSION 

In this section, we devote a short discussion to each of three 
important topics that relate to L/IMXB evolution. Binary mil- 
lisecond pulsars are discussed, with emphasis on the orbital pe- 
riod distribution and the long-standing birthrate problem. We 
also address the possible importance of X-ray irradiation of 
the donor stars in L/IMXB. Finally, we briefly investigate the 



possibility of forming low-mass black holes in L/IMXBs via 
accretion-induced collapse. 

5.1. Binary Millisecond Radio Pulsars 
and the Birthrate Problem 

Binary millisecond radio pulsars (BMPs) are widely thought 
to be the evolutionary descendants of L/IMXBs (Alpar et al. 
1982; Joss & Rappaport 1983). Using our BPS code and li- 
brary of evolutionary sequences, we compute the orbital period 
distribution of BMPs in the following way. Most sequences 
are terminated after the He or HeCO white-dwarf core of the 
donor star is exposed (see Paper I). For the binaries that con- 
tract to P or b < 2 hr, the evolution ends not long after the period 
minimum is reached; as mentioned above, further evolution to 
longer periods is not followed. At the end of each sequence, 
we know the orbital period and mass of the secondary. Each se- 
quence i and the corresponding final orbital period are weighted 
by the formation efficiency Tu and the results are accumulated 
to generate a histogram. We have not attempted to include any 
estimates for the pulsar lifetime or detectability in our analysis. 

Two examples of the calculated orbital period distribution are 
shown in Fig. 3 (solid histogram), where y = 0.0, ij M b = 10, 
at = 200kms _1 , and Ace = 0.1 and 0.5. Systems in Fig. 3 
with P 01 b < 2hr are meant only to indicate the relative pro- 
portions of short-period and long-period BMPs produced in 
our simulations. Overlayed is the period distribution of the 
observed systems (dashed histogram; taken from Taam et al. 
2000). Both distributions are normalized to unit area. Clearly, 
neither parameter set adequately reproduces the observed dis- 
tribution, and the agreement is extremely poor for Ace = 0.1. 
Larger values of Ace are strongly favored, since L/IMXBs form 
over a much wider range of initial periods and donor masses 
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(see Fig. 1, top panels) than when Ace ~ 0.1, which results in 
a wider range of BMP orbital periods. The same conclusion 
was reached by Willems & Kolb (2002), who used simplified 
models of L/IMXB evolution in a focused population synthesis 
study of BMPs. However, Ace ~ 0.5 is problematic, since the 
number of luminous L/IMXBs at the current epoch is greatly 
overproduced, by factors of >100, relative to the observed 
number, as noted in the last section. This problem emerges 
again when one considers the birthrates of BMPs based on the 
observed sample. 

Lorimer (1995) and (Cordes & Chernoff 1997) each con- 
servatively estimate the total Galactic birthrate of BMPs to be 
~10~ 6 yr _1 . Their likelihood analyses considered ^20 BMPs 
and included known pulsar selection effects, a model for the 
spatial distribution of BMPs, and estimated distance errors. 
Lorimer (1995) notes that several uncertainties may increase 
the birthrate by a factor of 10, to ^10" 5 yr _1 . In addition, the 
pulsar spin-down ages used in these studies may systematically 
overestimate the true ages of the systems (e.g., Lorimer 1995; 
Hansen & Phinney 1998), leading to larger actual birthrates. 
Given the ^100 observed LMXBs in the Galaxy and a typi- 
cal observable X-ray lifetime of (r x ) ~ 10 9 yr, the the semi- 
empirical birthrate of LMXBs is ~10~ 7 yr _1 , which is 10-100 
times lower than the BMP birthrate. Kulkarni & Narayan 
(1988) were the first to point out this potentially important dis- 
crepancy. 

Recall from § 3 that the birthrates of L/IMXBs lie in 
the range ^10~ 6 -10~ 4 yr _1 , which covers the range of semi- 
empirical BMP birthrates. The largest L/IMXB birthrates, 
which may be required if the BMP birthrate is ~40~ 5 yr _1 , re- 
sult from Ace ~ 0.3-0.5, which, in turn, yields far too many 
luminous X-ray sources at the current epoch. Although the ac- 
tive X-ray lifetime is not relevant for our theoretical L/IMXB 
birthrate calculations, this is an important quantity to calculate. 

We define the mean luminous X-ray lifetime from our simu- 
lations to be 

(rx) = ^^Af,-/^^, (12) 

where Af, is the total time spent by sequence i in the interval 
M a > 10~ 10 MQyr _1 , and the sums are over all sequences with 
Md < 4 selected in the population synthesis calculation. We 
have not attempted to correct (r x ) for systems that are transient 
according to the disk instability model, since the transient duty 
cycle is unknown. From Table 1 we see that (rx) is consistently 
^5 x 10 8 yr. This justifies the past use of (rx) ~ 10 9 yr as a typ- 
ical LMXB lifetime. However, it is essentially this long X-ray 
lifetime that yields such large numbers of luminous sources at 
the current epoch in our model calculations. We now discuss a 
possible resolution to this serious conflict. 

5.2. Irradiation-Induced Mass-Transfer Cycles 

Our binary evolution calculations do not account for X-ray 
irradiation effects on the secondary, which can dramatically 
change the evolution of the system, either by driving winds 
(Ruderman et al. 1989) or the expansion of the star (Podsiad- 
lowski 1991). Cyclic mass transfer may result from the irradia- 
tion, characterized by short episodes of enhanced mass transfer 
and long detached phases (e.g., Hameury et al. 1993; Harpaz 
& Rappaport 1994). Such cycles may then significantly reduce 
the total X-ray active lifetime, and thus resolve the L/IMXB 
overproduction problem and the discrepancy between the semi- 
empirical birthrates of BMPs and LMXBs. 



i 1 1 1 1 1 1 1 1 1 1 1 r 




cr k = 200 km s" 1 
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log Accretion Rate (M yr -1 ) 

FIG. 4. — Cumulative distribution of accretion rates at the current epoch 
for the same systems used to generate Fig. 2 (A C e = 0.5). Here we include 
both persistent and transient systems. The thick line corresponds to our ad hoc 
inclusion of irradiation effects. 

X-ray irradiation affects stars of mass <1.5M by ionizing 
the hydrogen at the base of the irradiated surface layer and dis- 
rupting the surface convection zone. By changing the surface 
boundary condition in this way, the star would like to expand to 
a new thermal-equilibrium radius (Podsiadlowski 1991), with 
more dramatic expansion as the total stellar mass is decreased 
and the fractional mass in the surface convection zone is in- 
creased. Intermediate-mass stars in IMXBs, which are initially 
fully radiative, may be strongly affected by irradiation once the 
mass is reduced to <1-1.5M Q and a surface convection zone 
appears, after the early, rapid phase of thermal-timescale mass 
transfer. 

An irradiation temperature of >10 4 K — a characteristic ion- 
ization temperature for hydrogen — gives a critical X-ray flux 
of S c ~ lO'^lO^ergss^cm" 2 , above which irradiation is 
important. Letting S c .u be the critical flux in units of 
10 1 1 ergs s" 1 cm" 2 , and M a _g be the accretion rate onto the NS in 
units of 10" 8 M Q yr _1 , we estimate the maximum orbital period 
for which irradiation is important: 

P orb ~70day(^) 3/ \ (13) 

where e < 1 is a factor that takes into account the geometry of 
the accretion disk and star, albedo of the star, and fraction of 
X-rays that penetrate below the stellar photosphere (Hameury 
et al. 1993). The value of e is uncertain, as is the effect of irra- 
diation on stellar and binary stellar evolution. However, it is at 
least quite plausible, and perhaps likely, that X-ray irradiation 
significantly alters the evolution of most L/IMXBs. 

When irradiation effects are relatively moderate, the results 
of Hameury et al. ( 1 993) indicate that irradiation-induced mass- 
transfer cycles do not typically change the secular evolution 
of P or b or Md- Under these circumstances, the inclusion of 
irradiation would not change the P or b and Md distributions in 
Fig. 2. However, the distribution of M a would change signifi- 
cantly. Suppose that the donors in all L/IMXBs are affected by 
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FIG. 3. — Orbital period distribution of binary millisecond radio pulsars. The solid and dashed histograms are for the theoretical and observed systems, respectively. 
Both distribution are normalized to unit area. 



irradiation in the same way, such that the mass-transfer rate is 
on average enhanced by a factor of f en > 1 during episodes of 
Roche-lobe overflow. We then expect the logM fl distribution in 
Fig. 2 to be shifted to the right by an amount log f en , but still cut 
off at the Eddington limit. Furthermore, the distribution is ev- 
erywhere decreased in height by an amount -log/ en , since the 
X-ray lifetime for each system is reduced by a factor f~ l . In 
order to alleviate the L/IMXB overproduction problem and the 
BMP birthrate problem, we perhaps need to have (rx) reduced 
by a factor of >100. This would require / en > 100, depending 
on the shape of the M a distribution calculated without irradia- 
tion. 

In Figs. 2 (Ace = 0.5) and 4 we have overlayed (thick, solid 
histogram) a distribution that may very roughly illustrate the 
net effects of irradiation for an illustrative enhancement factor 
°f /en = 30. To get this result, we simply multiplied the mass- 
transfer rate from the donor by / en and all times by in each 
of the selected evolutionary sequences, regardless of the instan- 
taneous values of Md or P 0I \,. We have chosen / en = 30 since 
(rx) and Nx are reduced by a factor of 10. The resulting M a dis- 
tribution (Fig. 2) is confined to the range M a > 10- 10 M Q yr- 1 , 
and we find that most of these systems are persistent. From 
the corresponding cumulative distribution shown in Fig. 4, we 
find that the ratio of the number of LMXBs with luminosities of 
10 36 -10 37 ergss _1 to the number with L x > 10 37 ergss _1 is ~0.6, 
smaller than the observed value (see § 4.4). However, irradia- 
tion effects at least produce the desired outcome of reducing 
this ratio from our earlier quoted theoretical value of <~5.3. Of 
course, our treatment of irradiation effects is greatly oversim- 
plified, and much more detailed calculations are requited. 

5.3. Low-Mass Black Holes 

It is generally assumed that the NS cannot accrete at rates 
exceeding the Eddington limit (M Edd ~ 10~ 8 M Q yr" 1 ). Material 
that is donated faster than M Edd , such as occurs during thermal- 
timescale mass transfer, will likely be ejected from the system, 



possibly as a radiatively driven wind from the accretion disk or 
in the form of relativistic jets. Evidence for both of these pro- 
cesses is seen in the X-ray binary SS 433 (Blundell et al. 2001), 
a system known to be in a phase of super-Eddington mass trans- 
fer. 

Our library of evolutionary sequences were computed with 
one reasonable, though heuristic, prescription for the mass cap- 
ture fraction (3: 



(b , \M 2 \<M Edd 
\bM Edd /\M 2 \, \M 2 \>M Eddl 



where b < 1 is a constant. This formula limits the accretion rate 
to be M a < bM^. We used a specific value of b = 0.5 in our 
calculations. 

The distribution of final NS masses (Fig. 5) was computed 
in the same way as the orbital-period distribution for BMPs. 
Masses up to ^2.5 M Q were reached in these simulations. Most 
modern NS equations of state give maximum masses of ^2- 
3M Q . Thus, even for b = 0.5, it is possible that a significant 
fraction of L/IMXBs ultimately contain a low-mass black hole. 
If we had used b = 1, our NS mass distribution would have 
broadened to include masses up to ^3.6M Q , and perhaps sev- 
eral tens of percent of NSs in L/IMXBs would collapse to black 
holes. No binary that contains a compact object of mass ^2- 
3M Q has yet been confirmed observationally. Such systems 
would not exhibit X-ray bursts and are unlikely to show the 
twin kHz quasi-periodic oscillations seen in LMXBs (van der 
Klis 2000). The inclusion of X-ray irradiation on the donor star, 
and the possible dramatic decrease in the X-ray lifetime, may 
generally yield final NS masses much closer to the initial mass 
of 1.4M Q . Perhaps this could explain why the measured NS 
masses in BMPs lie in a narrow range about ~ 1 .4M Q (Thorsett 
& Chakrabarty 1999). 
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FIG. 5. — Distribution of final NS masses for BPS parameters y = 0.0, 
Vmb = 1, (J* = 200 kms~', and Ace = 0.5. 

6. SUMMARY AND OUTLOOK 

Here we list the key points and results of this paper, in de- 
creasing order of importance. 

1. This is the first population synthesis study of L/IMXBs 
that incorporates detailed evolutionary calculations. With this 
addition, we are able to follow a population of L/IMXBs from 
(i) the incipient stage, (ii) to the current epoch, and finally (iii) 
to the remnant state when they presumably become BMPs. We 
are thus able to meaningfully compare our results with the sam- 
ple of observed LMXBs and BMPs. 

2. We have demonstrated that incipient IMXBs outnumber 
incipient LMXBs typically by a factor of >5 (see Table 1). 
Since IMXBs may evolve to resemble observed LMXBs, we 
claim that the majority of observed systems may have started 
their lives with intermediate-mass donor stars. 

3. We find that rather large values of Ace (~0.5) are re- 
quired in order (i) for the theoretical BMP P or b distribution to 
even remotely resemble the observed distribution, and (ii) to 
yield L/IMXB birthrates which are consistent with the semi- 
empirical BMP birthrates. However, we have discovered that 
such values of Ace lead to a dramatic overproduction of the 
number of luminous X-ray binaries in the Galaxy at the current 
epoch, by factors of - 100-1000. 

4. The overproduction problem and the discrepancy between 
semi-empirical BMP and LMXB birthrates may be resolved if 
the mean X-ray lifetime of L/IMXBs is reduced by a factor of 
> 100. Cyclic mass transfer, induced by the X-ray irradiation of 
the donor star, may have the desired outcome. It is plausible X- 
ray irradiation strongly affects the evolution of most L/IMXBs. 

5. Eddington-limited accretion onto the NS can lead to large 
NS masses of ^2-AM Q . It is then possible that a significant 
fraction of NSs in L/IMXBs collapse to low-mass black holes. 



To our knowledge, the only way to confirm the presence of a 
low-mass black hole in an observed LMXB is to measure the 
mass dynamically, which is extremely difficult in general. If the 
X-ray lifetimes of L/IMXBs are reduced substantially, by, e.g., 
the effects of X-ray irradiation, so too is the amount of mass 
that NSs may accrete. This may also explain why the measured 
NS masses in observed BMPs lie near <~1.4M Q . 

We conclude by listing a number of ways that our work may 
be extended and improved. 

1 . Since X-ray irradiation of the donor star may be an ex- 
tremely important component of L/IMXB evolution, it is im- 
portant to have a better quantitative understanding of this pro- 
cess. The problem is inherently three-dimensional, and must 
be treated as such in order to obtain meaningful quantitative re- 
sults. Important first steps in this regard have been made by 
Beer & Podsiadlowski (2002), but much work remains to be 
done. 

2. A detailed comparison of theoretical models of L/IMXBs 
formation and evolution with the observed population requires 
that we compute the spatial trajectories of the synthesized sys- 
tems in a realistic Galactic gravitational potential. It then be- 
comes possible to generate theoretical X-ray flux distributions 
as well as distributions in Galactic latitude and longitude. A 
treatment of observational selection effects, such as instrumen- 
tal flux limits and X-ray absorption, also becomes possible. 

3. Our library of L/IMXB evolutionary sequences will soon 
be incorporated into a sophisticated dynamical Monte Carlo 
code to study the evolution of globular clusters (Fregeau et al. 
2003, and references therein). The code now incorporates direct 
numerical integrations of single and binary-binary dynamical 
interactions, a mass spectrum of stars, and analytic treatments 
of single-star evolution. The inclusion of our library of X-ray 
binary calculations will make it possible to study directly the 
LMXB and BMP populations in globular clusters. 

4. An appropriately detailed treatment of the tidal evolu- 
tion of incipient L/IMXBs immediately after the SN should 
be included in future studies. Orbital circularization and spin 
synchronization of the secondary, coupled with orbital angu- 
lar momentum loss due to GR, as well the loss of spin angular 
momentum that results from MB should be considered. 

5. It would be extremely advantageous if L/IMXB evolution- 
ary calculations could be carried out at least 100 times faster 
than is possible at present. On a reasonably fast workstation, a 
typical computing time is currently ^20min, so that, realisti- 
cally, several days of computing time would be required to re- 
generate our current library of sequences. The development of 
an ultrafast Henyey-type stellar evolution code, along with in- 
creased processor speed over the next several years may make 
it possible to carry out 1000 L/IMXB evolutions in ^10hr of 
computing time. 
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